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SUMMARY 


The  results  of  a  series  of  testa  to  determine  the 
effect  of  internal  pressure  on  the  buckling  load  of  a 
thin  cylinder  under  an  applied  torque  indicated  that 
internal  pressure  raises  the  shear  buckling  stress. 

The  experimental  results  were  analyzed  with  the  aid  of 
previously  developed  theory  and  a  simple  interaction 
formula  was  derived. 


INTRODUCTION 


The  curved  metal  skin  of  a  modern  airplane  in 
flight  is  subject  to  stresses  that  may  cause  the  skin 
to  buckle,  and  proper  design  of  airplane  structures 
requires  a  knowledge  of  the  stress  conditions  under 
which  buckling  wili  occur.  The  ability  to  estimate  the 
buckling  point  under  combined  loading  conditions  is  of 
particular  importance. 

In  order  to  determine  the  effect  that  normal  pres¬ 
sure,  or  air  loads,  might  have  on  the  critical  stress 
for  curved  sheet,  two  preliminary  tests  were  made  (ref¬ 
ences  1  and  2) .  A  pronounced  increase  in  critical 
stress  with  increase  in  normal  pressure  was  found,  and 
the  subsequent  interest  shown  by  the  aircraft  industry 
in  this  subject  indicated  the  desirability  of  further 
study. 

No  well-established  theories  for  buckling  of  curved- 
sheet  panels  either  in  torsion  or  under  hydrostatic 
pressure  are  available;  however,  satisfactory  theories 
for  buckling  of  complete  cylinders  under  these  loadings 
have  been  advanced  (references  3  and  4) .  In  order  to 
effect  some  correlation  between  theory  and  experiment. 
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therefore,  an  investigation  was  made  of  the  influence  of 
•internal  pressure  on  the  critical  stresses  of  thin 
cylinders  in  torsion.  Experiments  were  conducted  to 
determine  the  critical  shear  stresses  for  four  cylinder 
lengths  at  a  number  of  different  internal  pressures. 

The  theories  and  the  experimental  data  were  used  in  con¬ 
junction  to  determine  an  interaction  formula  for  the 
buckling  of  cylinders  of  moderate  length  under  the  com¬ 
bined  action  of  torsion  and  internal  pressure. 


SYMBOLS 


E 

G 

Lo 

L 

t 


Young’s  modulus  of  elasticity,  psi 
shear  modulus  of  elasticity,  psi 
length  of  cylinder  without  rings,  inches 
length  of  cylinder  between  rings,  inches 
thickness  of  cylinder  wall,  inches 


d  diameter  of  cylinder,  inches 

r  radius  of  cylinder,  inches 


P 


t=0 


internal  pressure  of  cylinder,  considered 

positive  when  it  produces  tensi  le  stresses 
in  cylinder  walls,  psi 

critical  pressure  in  the  absence  of  torsion 
(negative  according  to  sign  convention 
adopted  for  p),  psi 


T 


shear  stress  in  cylinder  walls  due  to  applied 
torque,  psi 


Tcr  value  of  shear  stress  when  cylinder  is  at  the 

point  of  elastic  instability,  psi 

^Tcr)  critical  shear  stress  in  the  absence  of 

P“°  internal  pressure,  psi 

pressure  ratio 
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q 

e 

Ks 


7*7! 

V- 


shear-stress  -ratio  .. 


exponent  of  Rs  in  the  interaction  formula 

rotation  of  free  end1  of  cylinder,  radians 

coefficient  used  in  Lundquist's  empirical 
formula  (appendix  A) 

differences  in  strain  gage  readings 

Poisson's  ratio 


SPECIMENS  AND  TEST  PROCEDURE 


Diagrams  of  the  test-cylinder  construction  and 
ring  systems  used  are  given  in  figure  1.  The  cylinder 
was  made  from  0.032-inch  24S-T  aluminum- alloy  sheet 
closely  riveted  around  two  heavy  steel  rings,  one  at 
each  end.  The  sheet  was  Joined  along  an  element  of  the 
cylinder  with  a  butt  Joint  covered  with  a  single  strap 
on  the  outside.  Rings  made  of  24S-T  aluminum  alloy 
were  added  to  this  cylinder  and  divided  it  into  shorter 
cylinders,  the  lengths  of  which  were  equal  to  the  ring 
spacings . 

Figure  2  is  a  photograph  of  the  apparatus  used  to 
test  the  cylinder.  The  ends  of  the**  cylinder  were 
closed  by  heavy  flat  steel  plates  in  order  that  air 
could  be  maintained  under  pressure  inside  the  cylinder. 
The  weights  of  the  steel  plate  and  test  ring  at  the 
free  end  of  the  cylinder  were  neutralized  by  an  upward 
load  on  one  of  the  torque  arms.  Rotations  of  the  free 
end  of  the  cylinder  relative  to  the  floor  were  measured 
by  a  pair  of  dial  gages. 

Observed  buckling  loads  at  zero  and  at  low  internal 
pressure  were  determined  as  the  loads  at  which  there 
was  a  sudden  snap  of  the  cylinder  to  the  buckled  state. 
This  snap  action  was  in  many  cases  preceded  by  a  slow 
growth  of  visible  wrinkles  in  the  cylinder  walls  as  the 
load  increased  and  a  greatly  increased  rate  of  growth 
of  the  wrinkles  at  close  proximity  to  the  snap-buckling 
load.  With  increases  in  the  internal  pressure,  the 
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snap  decreased  in  violence  until  it  was  no  longer 
observable.  In  these  casea  the  buckling  load  was  esti¬ 
mated,  on  the  basis  of  visual  observations,  to  be  that 
load  at  which  the  rate  of  growth  of  wrinkles  with  load 
was  comparable  with  the  rate  of  growth  of  wrinkles  just 
prior  to  buckling  in  those  cases  in  which  snap  buckling 
did  occur. 

The  Southwell  method  of  determining  the  critical 
stress  of  specimens  with  initial  eccentricities  was 
assumed  to  be  applicable  to  cylinders  subjected  to 
torsion  and  was  used  as  a  check  on  the  visual  determi¬ 
nation,  Two  electrical  strain  gages  were  mounted  on 
opposite  sides  of  the  cylinder  wall  at  the  crest  of  a 
buckle,  the  location  of  which  was  determined  by  means  of 
a  preliminary  buckling  test.  The  difference  y  in 
readings  of  the  strain  gages  on  opposite  sides  of  the 
cylinder  wall  provided  a  numerical  measure  of  the  dis¬ 
tortion  as  the  wrinkles  grew.  The  Southwell  method  was 
applied  by  plotting  (y  -  y-^/^T  -  against  y  -  y-^ 

where  y^  and  were  arbitrarily  chosen  initial 

values  of  each  quantity  (reference  5) .  The  inverse 

slope  of  the  straight  line  formed  in  this  manner  is 

t  -  T,  where  t  is  the  desired  critical  stress, 
cr  1  cr 


DISCUSSION  OF  RESULTS 


In  figure  3  the  shear  stress  calculated  from  the 
external  load  on  th£  cylinder  is  plotted  against  the 
corresponding  rotation  of  the  free  end  for  various 
values  of  ring  spacing  and  internal  pressure.  The 
solid  lines  give  rotations  computed  by  the  formula 


T  2L0 

3  ~3~ 


in  which  G  was  assumed  to  be  3,970,000  psi. 

Prior  to  buckling,  the  stiffness  of  the  cylinder  in 
torsion  (slope  of  the  curves)  is  essentially  unaffected 
by  the  internal  pressure  or  by  the  ring  spacing.  The 
rotations  measured  were  slightly  greater  than  those 
predicted  by  the  formula.  This  result  may  be  attributed 
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to  the  method  of  measurement,  which  considered  as  rota¬ 
tions  many  other,  small  effects,  such  as  bending  of  the 
column  supporting  the  cylinder,  distortion  of  rivets  and 
rivet  holes,  and  taking  up  slack  in  the  bolt  holes  in 
the  end  plates. 

The.  experimental  buckling  stress  in  torsion  without 
pressure  is  compared  with  the  predictions  of  Donnell 
(reference  3)  and  Lundquist  (reference  6)  in  the  fol¬ 
lowing  table: 


L/d 

Critical  shear  stress. 

ksl 

Experimental 

Lundquist 

Donnell 

0.18 

5.8 

6.0 

7.9 

.36 

3.9 

4.2 

4.9 

.72 

3.0 

3.0 

3.3 

1.43 

2.3 

2.3 

2.3 

Lundquist' s  formula  was  derived  for  stress  at  failure. 
Except  for  very  short  cylinders,  however,  the  shear 
buckling  stress  and  the  stress  at  failure  are  essen¬ 
tially  the  same.  Donnell's  theoretical  curve  (solid 
curve  of  fig.  4)  and  Lundquist’ s  empirical  formula  are 
discussed  more  fully  in  appendix  A. 

In  figure  5  the  buckling  stresses  under  torsion 
are  plotted  against  the  internal  pressure  for  each 
cylinder;  this  figure  indicates  that  the  buckling 
stresses  increase  as  pressure  increases.  The  results 
were  essentially  the  same  whether  the  visual  or  the 
Southwell  method  of  determining  buckling  loads  was 
used,  although  the  Southwell  method  usually  gave 
slightly  higher  loads  than  the  visual. 

The  buckled  cylinder  could  be  returned  to  the 
unbuckled  state  either  by  decreasing  the  torque  or  by 
increasing  the  pressure.  The  loads  at  which  the  buckles 
disappeared  were  determined  visually  for  cylinders  1, 

2,  and  3a,  and  are  also  shown  in  figure  5.  These  loads 
were  not  recorded,  however,  for  the  remaining  cylinders 
because  of  a  large  scatter  in  the  readings.  This 
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scatter  Is  attributable  to  the  fact  that,  as  the  cylin¬ 
ders  became  shorter,  the  disappearance  of  the  buckles 
became  more  gradual  and  the  visual  selection  of  a 
"point  of  disappearance"  became  difficult.  In  the  case 
of  cylinder  2,  in  which  the  readings  are  relatively 
definite,  the  location  of  the  curve  indicating  disap¬ 
pearance  of  buckles  was  the  same  for  either-  an  increase 
of  pressure  or  a  reduction  of  torsion. 

The  method  used  to  find  the  Interaction  formula 
best  representing  the  experimental  data  is  given  in 
appendix  A.  The  analysis  leads  to  a  formula  of  the 
type 


R„q  +  Rp  =  i 

where  the  value  of  the  exponent  q  depends  upon  the 
assumptions  made  for  (pcr)T-o  an<*  (Tcr)p-o*  ^or 

cylinders  of  moderate  length,  which  according  to 
Donnell’s  analysis  satisfy  the  inequality 


the  exponent  q  is  equal  to  1#89  to  2  using  Tcr 
derived  from  Donnell’s  theoretical  curve,  the  exact 

value  depending  on  the  value  of  L^/td  concerned. 

Yfhen  Tcr  derived  from  Lundquist's  empirical  formula 

is  used,  q  =  2.17. 

In  figure  6,  curves  representing  the  exponents 
1.89,  2,  and  2.17  are  drawn  through  the  experimental 
data  of  figure  5  replotted  by  a  method  explained  In 
appendix  A.  The  three  curves  lie  close  together  and, 
so  far  as  fit  of  data  is  concerned,  little  basis  exists 
for  a  choice  among  them.  Simplicity  and  proximity  to 
the  average  value,  however,  recommend  use  of  q  =  2. 

The  equation  then  becomes 

R82  +  Rp  =  1 
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The  parameter  L2/td,  which  according  to  theory 
determines i  buckling  behavior,  may  be  varied  by  changing 
the  length,  thickness,  diameter,  or  any  combination  of 
dimensions.  For  this  reason,  the  restriction  of  the 
tests  to  one  thickness  and  one  diameter  probably  does 
not  constitute  a  significant  loss  in  generality. 

The  question  of  the  applicability  of  the  formula  to 
curved  panels  is  discussed  briefly  in  appendix  B. 


CONCLUSION 


The  critical  stress  of  a  cylinder  in  torsion 
increases  as  the  internal  pressure  Increases.  The  fol¬ 
lowing  interaction  formula  was  found  to  represent 
approximately  the  buckling  of  a  cylinder  of  moderate 
length  under  the  combined  effects  of  torsion  and  internal 
pressures 


R 


s 


2 


+  RP  = 


1 


where  Re  is  the  ratio  of  critical  shear  stress  with 

internal  pressure  to  critical  shear  stress  without 
internal  pressure  and  Rp  is  the  ratio  of  internal 

pressure  to  the  critical  pressure  in  the  absence  of 
torsion. 


Langley  Memorial  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va. 
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APPENDIX  A 

DERIVATION  OP  INTERACTION  FORMULA 


Introductory  Discussion 


Any  attempt  to  determine  which  of  the  conventional 
types  of  interaction  formulas  best  fits  the  experimental 
data  presented  in  this  paper  is  complicated  by  the 
question  of  which  values  to  use  for  /p  v  and 

\  °*VT=o 

(Tcr)  q*  Experimental  vaiuea  for  (Tcr)p=o  were 

obtained,  but  the  corresponding  experimental  values  of 
(Pcr^T_Q  were  not  investigated  because  of  the  destruc¬ 


tive  effects  on  the  cylinder  involved  in  such  a  test. 

A  simple  experimental  result  makes  possible  a  derivation 
of  an  interaction  formula  without  necessitating  assump¬ 
tions  concerning  stress  ratios  or  the  type  of  formula 
to  be  used. 


The  experimental  data  appear  to  indicate  that  the 
interaction  curves  for  the  four  cylinders  tested  are 
identical,  differing  only  in  position.  (See  fig.  5.) 
The  effect  of  a  change  of  length  appears  to  be  a  shift 
of  the  curve  parallel  to  the  pressure  axis.  In  fig¬ 
ure  6  this  effect  is  made  clearer  by  plotting  all  the 
data  to  a  common  p-intercept.  The  curve  for  i  =  1.43 


is  shifted  parallel  to  the  p-axis  a  distance  equal  to 
“(]?cr)T— o  (as  given  by  equation  (5)  which  follows)  in 


order  to  make  the  p-intercept  zero.  Each  of  the  other 
curves  is  then  shifted  the  proper  distance  for  best 
superposition.  The  relatively  slight  scatter  shows 
that  the  experimental  curves  are  nearly  superposable . 


If  it  is  assumed  that  the  interaction  curves  for 
different  lengths  of  cylinder  are  for  practical  pur¬ 
poses  identical,  a  computation  of  the  equation  for  the 
curves  from  theory  is  possible.  Among  the  points  on 
the  curve  of  figure  6  are  four,  indicated  by  modified 
symbols,  representing  the  special  case  p  =  0  for  the 
four  values  of  the  length.  Because  the  curve  of  figure 
represents  Tcr  plotted  against  P  -  (Pci.}t_q»  these  four 

points  represent  (TCr)  _0  Pi°tteci  against  "(pcr)r-o 
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for  var.iou.3  values  of  the  length.  The  formula  for  this 
curve  can  therefore '  he  obtained-  by  eliminating  the 
length  from  the  equations  expressing  (Tcr)  _0  fl^d 


(Pcr)T=0  ln  terms  of  L»  d* 


and  t. 


The  interaction 


curve  for  a  cylinder  of  given  dimensions  can  then  be 
obtained  by  shifting  the  origin  to  make  the  p-intercept 
equal  to  (Per)  -q  ** or  that  fcylinder.  The  interaction 


formula  may  then  be  derived  by  finding  simplified  expres¬ 
sions  for  (TCr)p=o  and  (pcr)T=o'  eliminating  L, 

and  transforming  to  give  the  final  formula. 


Simplified  Expressions  for  (Tcr)p=o 

Two  simplified  expressions  for  (Tcr)p=:o  are 

derived,  the  first  based  on  Donnell’s  curve  (reference  3) 
and  the  other  on  Lundquist’s  formula  (reference  6). 

Donnell’s  theoretical  result  for  the  shear  buckling 
stress  of  thin-walled  cylinders  with  simply  supported 
ends  subjected  to  torsion  is  given  by  the  solid  curve  of 
figure  4.  For  cylinders  of  moderate  length,  the  curve 
is  nearly  a  straight  line  given  by  the  equation 


(la) 


while  for  extremely  large  values  of  L^/td  a  better 
fit  is  given  by  the  equation  of  the  asymptote 


(lb) 


This  essentially  straight  portion  of  the  curve  is 
included  between  the  limits 


30a  <  (h  f  <  sf 


(2) 
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At  the  lower  limit  of  Inequality  (2),  equation  (la)  Is 
In  error  by  about  15  percent  and  equation  (lb),  by  about 
20  percent.  (For  the  upper  limit  see  reference  3.) 

Lundquist’s  formula  Is 


KSE 


-1.35 


where  values  of  Kg  are  given  In  reference  6  for  various 
values  of  L/r.  In  order  to  determine  a  formula  for  Ka, 
the  values  of  Ks  were  plotted  against  L/r  on  loga¬ 
rithmic  paper.  (See  fig.  7.)  When  L/r  >  0.5, 


K 


s 


1.27 


(I)'0'46 


At  L/r  =  0.35  (approximately  the  lower  limit  of 
Inequality  (2)  for  the  cylinder  tested) ,  the  error  In 
this  formula  is  about  7  percent.  Use  of  this  value  of 
Kg  leads  to  the  formula 


cr 


p=0 


=  2.36E 


B(a) 


— 0 . 1 . 35 


(3) 


Simplified  Expression  for 


(Pcr)T=o 


A  formula  based  on  work  by  von  Mises  and  developed 
at  the  David  Taylor  Model  Basin  (equation  (10)  of  refer¬ 
ence  7)  for  the  buckling  of  a  closed  cylinder  under 
hydrostatic  pressure  may  be  written 


If,  as  In  the  Inequality  (2) , 
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the  denominator  of  equation  (4)  remains  within  10  percent 
of  unity,  so  that,  to  this  degree  of  approximation. 


“(pcr)T=0 


2.60! 


*$) 


,2.5 


(5) 


The  negative  sign  appears  because  In  the  present  paper 
Internal  pressure  Is  considered  positive. 


The  Interaction  Formula 

For  the  sake  of  generality  and,  at  the  same  time, 
of  simplicity  In  the  derivation  of  the  interaction 
formula,  equations  (1)  and  (5)  may  be  written 


Tcr)p=0  =  K2Lnf2(d,t) 

(6) 

Pcr)T=0  =  K1Lmf1(d,t) 

(7) 

where  and  are  arbitrary  constants,  m  and  n 

are  arbitrary  exponents,  and  f-^  and  fQ  are  arbitrary 

functions  of  diameter  and  thickness.  Elimination  of  L 
gives 


(Pcr)T=i0 


Klfl 


(Tcr)p=0 

K2f2 


m/n 


This  equation  is  the  formula  for  the  curve  of  fig¬ 
ure  6,  which  Is  a  plot  of  Tcr  against  p  -  ( Per) t-o 

and  therefore,  when  p  =  0,  a  plot  of  (TCr)  _q  against 

“(Pcr)T-o’  order  to  find  the  interaction  formula, 

substitute  p  -  (Pcr)T=0  in  Place  of  “(pcr)T=0  and 

Tcr  in  place  of  (Tcr)p-Q  an<i  simplify  the  resulting 

expression.  Thus, 
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P  -  (Pcr)T=0  = 


m/n 


Dividing  through  by  "(Pcr)r-o  using  equations  (6) 

and  (7)  yields 


..  P„  +  l  = 

(Por)T=0 


m/n 


cr 


~,m/n 


(Tcr)p=Oj 


A  rearrangement  of  terms  gives 

-.m/n 


cr 


C1*)p=0j 


(Per)** 


=  1 


If  Donnell's  equation  (la)  is  used  with  equation  (5), 

5  =  CT53  =  i-89 

or  if  Donnell's  equation  (lb)  and  equation  (5)  are  used, 

m  _  1  =  p 

n  **  0.5 

Use  of  Lundquist's  equation  (3)  with  equation  (5)  yields 

a  =  -i*  =  2.W 

n  0.46 

For  purposes  of  simplification,  let 


Rs  = 


cr 


(Tcr)p=0 


Rp  = 


(PCI»)T=g 
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and 


q.  = 


Then,  the  final  formula  may  be  written 


ejib 
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APPENDIX  B 

APPLICATION  OP  INTERACTION  FORMULA  TO  CURVED  PANELS 


The  theoretical  analysis  of  the  buckling  of  a 
curved  panel  is  much  more  difficult  than  the  treatment 
of  complete  cylinders.  In  the  case  of  buckling  under 
torsion  alone,  the  boundary  conditions  required  for 
panels  greatly  complicate  the  theoretical  analysis.  It 
may  be  presumed,  however,  that  the  added  restraints 
make  the  critical  stress  higher  than  that  of  a  complete 
cylinder  of  the  same  length,  thickness,  and  curvature. 

The  relation  between  the  pressure  and  the  direct  stresses 
produced  In  the  sheet  Is  also  much  more  complicated  than 
in  the  case  of  a  complete  cylinder.  The  axial  stresses 
depend  on  the  area  of  the  end  bulkhead  and  on  the  areas 
of  cross  section  of  sheet  and  3oar  elements.  It  is 
difficult  to  draw  conclusions  regarding  the  relation  of 
pressure  to  circumferential  stresses.  It  is  therefore 
not  to  be  expected  that  an  interaction  formula  established 
for  complete  cylinders  will  necessarily  apply  to  curved 
sheet  also. 

In  order  to  determine  whether  the  interaction 
formula  derived  in  appendix  A  may  be  used  to  gain  a  rough 
idea  of  the  strengthening  effect  of  internal  pressure  in 
the  case  of  curved  panels,  the  formula  is  applied  to  the 
test  points  of  the  two  specimens  of  reference  2.  Figure 
8  is  a  sketch  of  the  two  specimens  tested.  In  fig¬ 
ure  9  the  experimental  shear  buckling  stresses  are 
plotted  against  the  internal  pressure. 

Curves  representing  the  interaction  formula  are 
also  plotted  in  figure  9«  In  order  to  apply  the  formula, 
the  experimentally  determined  values  of  (Tcr)p=£  were 

used;  (Pcr)T-0  assUTr-e(5  to  be  the  same  as  for  a  com¬ 
plete  cylinder  of  identical  length,  thickness,  and 
curvature,  and  is  found  by  applying  equation  (it.)  of 
appendix  A.  Fair  agreement  was  obtained  between  the 
curves  and  the  test  points  and  the  agreement  might  be 
considerably  improved  by  a  more  precise  method  of 
determining  (PCr)  =0* 
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Figure  2.-  Test  setup. 
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Figure  3.-  Shear  stress -rotation  curves  for  cylinders 
with  various  ring  spacings  and  internal  pressures. 


Figure  4.-CriTica(  shear  stress  for  thin-wailed  cylinders. 
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Figure  5-  Effect  of  Internal  pressure  on  critical  shear  stress  for 
thin -walled  circular  cylinders. 
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Figure  7-  Value  of  Ks  used  in  empirical  formula  for 
critical  shear  stress  of  thin -walled  cylinders. 
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Figure  9.”  Interaction  formula  applied  to  curved  panel  specimens 
of  reference  2. 
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ABSTRACT: 


Tests  Indicate  that  Internal  pressure  raises  the  shear  buckling  stress.  Formula  was 
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internal  pressure  and  the  ratio  of  internal  pressure  to  the  critical  prejsure  In  the 
absence  of  torsion. 
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